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t:;^- ■ Color confinement can be understood by the dual Higgs theory, where monopole con- 

^ ' densation leads to the exclusion of the electric flux from the QCD vacuum. We study the 
role of the monopole for color confinement by investigating the monopole current system. 
When the self-energy of the monopole current is small enough, long and complicated 
monopole world-lines appear, which is a signal of monopole condensation. In the dense 
monopole system, the Wilson loop obeys the area-law, and the string tension and the 
monopole density have similar behavior as the function of the self-energy, which seems 
. that monopole condensation leads to color confinement. On the long-distance physics, 
^ ! the monopole current system almost reproduces essential features of confinement proper- 
^ \ ties in lattice QCD. In the short- distance physics, however, the monopole-current theory 
would become nonlocal and complicated due to the monopole size effect. This monopole 
size would provide a critical scale of QCD in terms of the dual Higgs mechanism. 
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1. QCD-monopoIe for Quark Confinement 



Quantum Chromodynamics (QCD) is the fundamental theory of the strong interaction 
I and describes not only hadrons but also their underlying structure in terms of quarks and 
gluons. In the high energy region, the QCD gauge coupling becomes small, and therefore 
the system can be described by the perturbative QCD, where quarks and gluons are 
good degrees of freedom. On the other hand, in the low energy region, there are various 
nontrivial phenomena like color confinement, dynamical chiral-symmetry breaking and 
topological defects due to the strong coupling nature of QCD. However, the essential 
degrees of freedom are unclear in the nonperturbative QCD, and their finding is strongly 
desired for understanding hadron physics. 

In 1974, Nambu proposed an interesting idea that quark confinement can be interpreted 
using the dual version of the superconductivity 0]. In this picture, condensation of the 
color magnetic monopole is the key concept and leads to squeezing the color-electric flux 
between quarks through the dual Meissner effect. In 1981, 't Hooft showed that QCD 
is reduced into an abelian gauge field theory with magnetic monopoles after the abelian 
gauge fixingp. In this gauge, the only abelian gauge fields with monopoles would be 
essential for describing non-perturbative phenomena in the low energy region|0-|^. 

Recent lattice QCD simulations show strong evidence on this dual Higgs picture for 
the nonperturbative QCD in the Maximally Abelian (MA) gauge[^-[^. The MA gauge 
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is considered as the best abehan gauge, and is reahzed by the gauge transformation so 
as to maximize the diagonal components of the gauge fields (or the link variables in 
the lattice formalism). In this gauge, information of the gauge configuration is con- 
centrated into the diagonal components. According to the nontrivial homotopy class, 
7r2(SU(A/'J/U(l)^=-i) = Z^^"^0, the QCD-monopole appears from the hedgehog-like 
configuration, which accompanies large fluctuation of the gauge field in the QCD vac- 
uum. Since the classical monopole mass is in inverse proportion to the gauge coupling 
constant, monopoles can be easily excited in the low energy region with the large gauge 
coupling. Therefore, QCD-monopoles are expected to be the relevant degrees of freedom 
for the description of the non-perturbative QCD. However, such a picture for the confine- 
ment mechanism in terms of monopole condensation can be applied only in the infrared 
region. In the ultraviolet region, the large monopole mass leads to suppression of the 
monopole excitation, and the QCD system seems to be trivial in terms of topology. 

In the dual Higgs phase, only short-range interactions remain among monopoles, and 
therefore the monopole current dynamics is subjected to a simple action. Thus, the 
nonperturbative QCD would be described only by the monopole degrees of freedom. In 
this paper, we study the monopole current dynamics and color confinement from the 
idealized monopole-current action. 

2. Monopole Current Dynamics and Kosterlitz-Thouless Type Transition 

We study the monopole current system using the idealized monopole action, and in- 
vestigate the role of the monopole currents for color confinement. In general, the current 
action includes nonlocal interaction terms. In the confinement phase, however, the nonlo- 
cal part of the interaction can be neglected, because the currents are expected to interact 
each other only in the short distance due to the screening effect by the dual Higgs mech- 
anism|^. With the current conservation, the partition functional is written as 

Z = y Dfc^exp{— a j d'^xtTk^{x)}S{d^kfi), (1) 

where fc^(x) = k^{x) -t^ /2 and a are the monopole current and energy for the unit current 
length, respectively. In order to perform the path integral (|l]), we put the system on the 
4-dimensional lattice with the lattice spacing a. In the lattice formalism, the monopole 
currents are defined on the dual lattice, = e/(47r) ■ a^k^{s). The partition function 

is given as Z = X^fciat exp{-a^''* J2s where a^""* = a/2 - (47r/ea)^. 

The lattice QCD simulation shows that one long monopole loop and many short loops 
appear in the confinement phase. The only long loop becomes important for the properties 
of the vacuum, while many short loops are originated from the fluctuation in the ultravi- 
olet region and therefore can be neglected. In this system, the partition function can be 
approximated as the single monopole loop ensemble with the length L, Z = J2l p{L)e~'^^, 
where L and p{L) are length of the monopole loop and its configuration number, respec- 
tively. The monopole current with length L is regarded as the L step self-avoiding random 
walk, where 2d — 1 = 7 direction is possible in each step. Therefore, p{L) is roughly es- 
timated as {2d — 1)^ = 7^. Using this partition function, the expectation value of the 
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monopole current length is found to be 

\ oo if a < ln(2(i - 1). ^ ' 

When the energy a is larger than the entropy ln(2(i — 1), the monopole loop length is 
finite. However, when a is smaller than the entropy, the monopole loop length becomes 
infinite, which is a signal of monopole condensation in the current representation. 

In performing the simulation, we construct the monopole current as the sum of pla- 
quettes of the monopole current because of the current conservation condition. First, we 
prepare a random monopole current system (hot start) and no monopole current system 
(cold start). Then, we update the link of the monopole current using the Metropolis 
method. 

We generate the monopole current system on 8^ lattices. Fig.l shows the monopole 
current vacuum for the typical cases, (a = 1.5, 1.8, 1.9) at a fixed time. Corresponding to 
each case, we show in Fig. 2 the histograms of the monopole loop length [p!0[] of one cluster 
with 40 configurations. For large a, only small loops appear and their density is small. 
On the other hand, for small a, there appear one long monopole loop and some short 
loops, and the monopole currents are complicated and dense. Between these two regions, 
there is critical value, Uc — 1.8, where monopole loops with various length appear. The 

monopole density pm = Z)s,/i |^Ai('5)| and the clustering parameter rj = ^-^^^ are 

shown in Fig.3. and Fig. 4, respectively. Here, Lj is the loop length of the z-th monopole 
cluster. The property of the clustering parameter is that it goes to unity when many loops 
cluster and unite to one long loop, and it approaches to zero when loops never cluster each 
other. When a becomes small, the monopole density becomes large gradually at a ~ etc- 
However, the clustering parameter is drastically changed at ac- Quantitatively, the critical 
value of monopole current energy, — 1.8, is close to ln7 ^ 1.95, which is corresponding 
to the entropy for the 1 step self-avoiding random walk. This vacuum can be roughly 
approximated with the vacuum made of the longest monopole loop alone in terms of the 
mechanism of the transition. Such a transition is quite similar to the Kosterlitz-Thouless 
type transition in 2-dimensional superconductors, where vortex condensation may occur. 

3. Dual Field Formalism and Role of Monopoles for Confinement 

In this section, we study how these monopole currents contribute to the color con- 
finement properties. Quark confinement is characterized by the linear inter-quark po- 
tential, which can be obtained from the area-law behavior of the Wilson loop, (H^) = 
{Pexp{ie § A^dx^)) . Therefore, it is desired to extract the gauge variable from the 
monopole current A;^. We now derive the abelian gauge variable in stead of the non- 
abelian gauge field, because the lattice QCD results show the color confinement phenom- 
ena can be discussed only with abelian part to some extent. However, in the presence 
of the magnetic monopoles, the ordinary abelian gauge field A^{x) inevitably includes 
the singularity as the Dirac string, which leads to some difficulties in the field theoretical 
treatment. Instead, the dual field formalism is much useful to describe the monopole cur- 
rent system, because the dual gauge field B^{x) can be introduced without the singularity 
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for such a system with /c^ 7^ and = 0. This is the dual version of the ordinary gauge 
theory with A^{x) for the QED system with 7^ and fc^ = 0. 

Let us consider the derivation of the dual gauge field B^{x) from the monopole current 
fc^(x). Using the relation d^,*F^y = with *F^^ = \e^yapFap, *F^,^ is given by the 2-form 
of the dual gauge field as *F^i, = d^By — dyB^, so that one finds d"^ By — dy{d^B^) = k^. 
By taking the dual Landau gauge, d^B^ = 0, this field equation becomes quite a simple 
form, d^Bfj^ = k^. Therefore, starting from the monopole current configuration kn{x), we 
derive the dual gauge field B^, 

BAx) = B-%(x) = - jij / d'y-^,. (3) 

Using the dual gauge field, the Wilson loop is obtained as (W) = {exp{ie § Afj_dx^}) = 
{exp{iele^^ap I da^^*Fap)) = {exp{iele^^af3 1 daf,„{daBf3 -dpB^)}). 

Now, we apply this dual field formalism to the monopole current system discussed 
in Section 2. Since the monopole current /cjf*(s) is generated on the lattice with the 
mesh a, the continuous dual field B^{x) is derived from k^{x) = An/e ■ — s) 

using Eq.(|^), in principle. In estimating the integral in Eq.(|^) numerically, we use a 
fine lattice with a small mesh c. To extract B^{x) correctly, the mesh c is to he taken 
enough small. However, too fine mesh is not necessary because the original current A;Jf*(s) 
includes the error in the order of a. Numerical analyses show that the use of c ~ a is 
too crude for the correct estimation of the integral in Eq.(|^). Instead, the calculation 
with c < a/2 is good enough for the estimation of B^{x), and hence we take c = a/2 
hereafter. On lattices, the dual gauge field 9'^^^^ = aeB^/2 is defined in the dual Landau 
gauge, 9^^^\s + fi) = 27iJ2s' A(s — s')A;|f*(s'), using the lattice Coulomb propagator A(s) = 
{dnd'^)~^, where d' denotes the backward derivative. The dual version of the abelian field 
strength Of^'^''^ = ea^*F^^/2 is defined by = ^^C^^'ls) - ^^^''Hs)- 

The expectation value of the Wilson loop (W) is shown in Fig. 5. The Wilson loop 
exhibits the area-law and the linear confinement potential: In(iy) decreases linearly with 
the inter-quark distance, where its slope corresponds to the string tension. Quantitatively, 
the string tension is measured by the Creutz ratio, and we show in Fig. 6 x(3, 3) as a 
typical example. For a < ac, the string tension decreases as a increases, and its behavior 
is similar to the lattice QCD result by setting a = 0.8/?. For a > Uc, the string tension 
becomes almost vanishes. It is worth mentioning that the behavior of string tension or 
the Creutz ratio becomes similarly to the monopole density shown in Fig. 3. Thus, the 
confinement force is controlled by the monopole density, and therefore monopoles can be 
regarded as the essential degrees of freedom for color confinement. 



4. Monopole Size and Critical Scale in QCD 

In the final section, we compare the monopole current system with abelian projected 
QCD I^J^ especially in terms of the coupling correspondence and the size of QCD- 
monopoles. 

To begin with, let us consider one magnetic monopole with an intrinsic radius R in 
the multi-monopole system. In the static frame of the monopole, it creates a spherical 
magnetic field, H(r) = = ior r > R and H(r) = = for r < R, 
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where the QCD running gauge couphng e(r) is used to include the vacuum polarization 
effect. Here, the effective magnetic-charge distribution is assumed to be constant inside 
the monopole for the simplicity. 

Now, we consider the lattice formalism with a large mesh a > R. The electromagnetic 
energy observed on the lattice around a monopole is roughly estimated as 



which is largely changed depending on the lattice mesh a. This simple estimation neglects 
the possible reduction of g{r) in the infrared region due to the asymptotic freedom of 
QCD. The screening effect of the magnetic field by other monopoles also reduces g{r) 
effectively in a dense monopole system. However, M{a) is modified by at most factor 2 
(M(a) = ^2 (^)q ) even for the screening case as g{r) = g{a) ■ 6{2a — r). Then, even in the 
multi-monopole system, M(a) would provide an approximate value for the electromagnetic 
energy on lattices created by one monopole, and we call M(a) as 'lattice monopole mass'. 

For the large mesh a > R, the monopole contribution to the lattice action reads 5* = 
M{a)a ■ L, where L denotes the length of the monopole current measured in the lattice 
unit a. Therefore, M{a) is closely related to the monopole-current coupling a^^^ and 
(3 = 2Nc/e'^, a^^^ ~ M(a)a ~ = | ■ /5su(2)- For the above screening case, this relation 
becomes a^^^ ~ M{a)a ~ = f ■ Psv{2), which is consistent with the numerical result, 
a = 0.8/?, discussed in Section 3. Here, as long as the mesh is large as a > R, the lattice 
monopole action would not need modification by the monopole size effect, and the current 
coupling a^^^ is proportional to (3. Quantum mechanically, there is the energy fluctuation 
about at the scale a, and therefore monopole excitation occurs very often at the long- 
distance scale satisfying M{a) ~ a'^. Thus, M{a)a ~ a^*^* is the control parameter for 
the monopole excitation at the scale a{> R), and we can obtain the quantitative criterion 
for 'monopole condensation' as M{a)a ~ ln{2d — 1) from the analysis using the current 
dynamics in Section 2. 

Second, we discuss the ultraviolet region with a < R. In the current dynamics, there 
exists a critical coupling ac — ln(2(i — 1) as shown in Section 2 and 3. Above ac, the lattice 
current action provides no monopole condensation and no confinement, while /3 — 00 can 
be taken in the original QCD keeping the confinement property shown in Fig. 6. Such a 
discrepancy between (3 and a can be naturally interpreted by introducing the monopole 
size effect. Obviously, the monopole-current theory should be drastically changed in the 
ultraviolet region a.s a < R, if the QCD-monopole has its peculiar size R. 

Here, let us reconsider the relation between a and a^*^* in the lattice current action. 
Similarly in the lattice QCD, the action has no definite scale except for the lattice mesh 
a(> R), and therefore the scale unit a would be determined so as to reproduce a suitable 
dimensional variable, e.g. the string tension a ~ IGeV/fm, in the monopole current 
dynamics. For instance, a is determined as a function of a'*^* using the Creutz ratio 
X — (JO? in Fig. 6. Therefore, a should reach R before realizing a^^* ajf*, and the 
framework of the current theory is to be largely modified due to the monopole size effect 
for a < R. 

Finally, we consider the fine lattice with a small mesh a < R, and study the monopole 
size effect to the current action (0). In this case, the extended structure of the QCD- 




(4) 



6 



monopole is observed on the lattice, and the monopole creates the electromagnetic energy 



Mia) . ,'.-Hirr - + drr^f^ = l^F = (5) 

on the lattice. Then, the 'lattice monopole mass' M(a) is determined mainly by the size 
R, and is almost a-independent. Accordingly, the lattice action should be changed as 

S = M{R)a k,{s)k,{s')-9(--\s-s'\)6s^,s'^ (6) 

because of the self-avoidness originating from the monopole size R. Thus, the monopole- 
current theory becomes nonlocal in the ultraviolet region a < R. 

In conclusion, the QCD-monopole size R provides a critical scale for the description 
of QCD in terms of the dual Higgs mechanism. In the infrared region as a > R, QCD 
can be approximated as a local monopole-current action^], and the QCD vacuum can 
be regarded as the dual superconductor. On the other hand, in the ultraviolet region as 
a < R, the monopole theory becomes nonlocal and complicated due to the monopole size 
effect, and the perturbative QCD would be applicable instead. 
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Figure Captions 

Fig.l The monopole currents for typical cases, (a) a = 1.9, (b) a = 1.8, (c) a = 1.5. 
Fig. 2 The histograms of the monopole current length of each cluster with 40 configura- 
tions for the same cases with Fig.l. 

Fig. 3 The monopole density pm as the function of a. The monopole density becomes 
small for large a and almost vanishes for a > Oc- 

Fig. 4 The clustering parameter t] as the function of a. At Oc, f] is drastically changed. 
Fig. 5 The expectation value of the Wilson loop {W{I x J)) on fine mesh c = a/2 for 
a = 1.7, 1.8, 1.9. The area-law behavior indicates the linear interquark potential. 
Fig. 6 The Creutz ratio as the function of a in the multi-monopole system. The dotted 
line denotes the Creutz ratio as the function of a = 0.8/3 in the lattice QCD. 
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